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abstract 

We investigate the critical dissipation of the double-well quantum mechanics. We adopt two-state 
approximation to define effective Ising models and apply the block decimation renormalization group 
and the finite range scaling method recently proposed for the long range Ising model. We briefly 
report the numerical results of the critical dissipation for various model parameters. 



1 Introduction 

The phase transition here is the quantum-to-classical phase transition. Consider a quantum 
system of the double well potential. Due to the tunneling, the system is oscillating and sym- 
metric, which we call 'quantum'. With dissipative effects, the tunneling is suppressed, and the 
ground state becomes a localized state to break the original Parity symmetry, which we call 
'classical'. Then, there must be a 'phase transition' between these two states. Our subject is 
to evaluate the critical dissipation for this phase transition using a new renormalization group 
method and a new scaling ansatz according to the range of interactions. 

The above picture about the phase transition has been common now by many researches 
with various methods [TJ [TBI El IB] • This type of transition is expected to be observed as deco- 
herence effects in the recent quantum mechanical few body experiments[10j or future quantum 
devices. Quantitative treatment of decoherence processes is essential in planing realistic quan- 
tum devices, but it must be a difficult subject since it need to bridge between the quantum 
and classical notions. Certainly the double-well quantum mechanics with dissipative effects is 
a first step model which we must fully understand quantitatively. 

The actual value of the criticality for the model as a function of the parameters of the double- 
well potential and the dissipative interactions, are not well established. Also, the double-well 
system has been frequently compared to the one- dimensional Ising model especially for the deep 
well case. However, the detailed and quantitative relations between the double-well quantum 
mechanics and the Ising model with general type of interactions have not been clear yet. 

In this article, we develop new methods of handling the infinite range interactions, that is, 
the block decimation renormalization group and the finite range scaling. These two methods are 
new way of evaluating the infinite range system with recourse to the intermediate finite range 
systems. The block decimation renormalization group is a real space renormalization group 
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method to solve finite range system exactly. Then the finite range scaling ansatz is applied 
to the susceptibility calculated exactly for finite range systems, and we evaluate its scaling 
exponent. The finite range scaling exponent finally tells us the value of the critical coupling 
constant through the pole position of the zeta function. These new methods may shed light 
on other infinite range systems even in higher dimensional space, and may open new analyzing 
power in the renormalization group world. 

First of all we should recall that it is non-trivial to include dissipative effects in the quantum 
mechanics, since there is no simple Hamiltonian to realize friction proportional to the velocity. 
We take the Caldeira-Leggett (CL) modelpLj, where dissipation comes out of the microscopic 
origin. The model consists of a target system and environmental degrees of freedom of infinitely 
many harmonic oscillators, 



S[ q, {x a }} =Jdt{ l -Mq 2 - V (q) + £ 



^m a x 2 a - ]^m a u 2 a x 2 a - qC a x c 



(1) 



The target variable q and each environmental variable x a are coupled linearly with coupling 
constant C a . Due to these couplings, the energy of the target system is transfered to the 
environmental oscillators, and it exhibits as dissipative effects. 

To describe the effective target dynamics, we path integrate out the environmental degrees 
of freedom. Then the non-local interactions of the target variable emerge as, 
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where 77 is an effective coupling constant to represent the strength of dissipation and parameter 
p denotes the damping rate of the non-local interactions. Of course these parametrization is 
a simplified expression of non-local interactions which are determined by micro parameters 
m a ,u a and C a . It is known that with quadratic damping (p = 2) interactions, the classical 

dq 

equation of motion of q(t) suffers effectively a friction proportional to the velocity — , and thus 

at 

this case is called the Ohmic case. 

In the Euclidean path integral formalism, the quantum mechanical system is equivalent to 
a one- dimensional statistical system. If we approximate the system using the minimum degrees 
of freedom, that is, two states par site (after discretization in time with a finite slice), it is 
transformed into the Ising spin chain with long range interactions. 

Long range Ising models with the above type of simple interactions, 



~ — ^2^2 K n a k a k+n = X! ~ a k^k+n , er fc = ±1 



nP 



(3) 



have its own long history pioneered by Refs.jU El [3], and it has been understood that with 
strong enough long range interactions there occurs the spontaneous magnetization. The known 
facts [SI El E] about the phase transition are that the critical coupling constant rj c is finite 
only for 1 < p < 2, that is, r] c = for p < 1 and rj c = 00 for p > 2. Particularly the p = 2 
case, which is called as "Ohmic" because it corresponds to the velocity proportional friction 
in the classical mechanical situation, is a very special boundary where it is proved that the 
Kosterlitz-Thouless transition occurs [9]. 

It is difficult, however, to evaluate the value itself of rj c . In the previous article we successfully 
evaluated the critical coupling constant of this long range Ising model, where we developed a new 
method of investigating the infinite range system, Finite Range Scaling (FRS) hypothesis [11] . 
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Figure 1: BDRG procedure 



In this article we apply this method to the double-well quantum mechanics. Here we briefly 
report the prototype calculation using two-state approximation of the system. Even with two- 
state approximation, the effective Ising spin interactions are not a straightforward mapping of 
the original non-local interactions. We will confirm that FRS method does work to evaluate 
the critical dissipation and our results are consistent with those obtained by other methods. 
The full calculation without state reduction will be reported elsewhere. 



2 Finite Range Scaling method in the long range Ising 
model 

For the long range Ising models defined in Eq.flS]), we successfully evaluated the critical dissipa- 
tive interactions to bring about the spontaneous magnetization. We developed a new method 
of the Block Decimation Renormalization Group (BDRG) and the Finite Range Scaling (FRS) 
ansatz [TTJ. Our method utilizes simple calculation with very small size of the computational 
resources. There are some large scale Monte Carlo studies [T2| [T3| [H] for the limited parameter 
space, and their results are perfectly consistent with ours, which indicates our ansatz of FRS 
is just on the right way of accessing properties of infinite range interactions. 

We explain the BDRG method in Fig.l in case of range n = 3. The DRG (Decimation 
Renormalization Group) is a block spin transformation method formulated by Wilson to cal- 
culate the partition function approximately [15]. In one dimensional nearest-neighbor Ising 
system, it works without any approximation and we can get the exact partition function very 
easily. In the case that there are non-nearest-neighbor interactions, the original DRG does not 
work. We developed an extended version, the Block DRG, which may accommodate long range 
interactions. 

We consider a system with long range interactions with coupling constants K m which couples 
m-separated two spins. We set the maximal finite range of interactions to be n, and investigate 
the infinite range limit by increasing n, one by one. First, we divide spins into blocks of size n. 
Then there are interactions only between nearest neighbor blocks, that is, looking at the model 
'block-wisely', it is nothing but the nearest-neighbor model where one block plays a site of 2 n 
states. 

Next, all the spins in the middle block in Fig.l are decimated, namely, are integrated for 
all possible 2 n states. After decimation we obtain new inter-block interactions which still keeps 
nearest-neighbor-block property. Now the renormalization transformation is defined and can 
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be repeated without approximation. 

We represent inter-block interactions by the transfer matrix T of 2 n x 2 n . Then the renor- 
malization transformation is nothing but to make a product of T matrices, 

y(fc+l) _ rp(k) . j.(fe) ^ ^ 

where matrix is inter-block interactions between /c-th renormalized neighboring blocks. 

By BDRG we calculate the susceptibility of the system xi n ) w ith a finite range interactions 
exactly. We impose the extend field h as, 



^2 K n a k a k+n + ha k 



k 

Then the susceptibility of the system is given by, 

1 d 2 (lnTrTW) 



(5) 



x(n) — hm 



fc->oo n ■ 2 k dh 2 



(6) 

h=0 



We increase the range n by 1 to see the change of the susceptibility and define the corre- 
sponding scaling exponent j3 as follows, 

/ 1 \ 0(p,v) 

logx(n) -logx(n-l) oc (^-J . (7) 

We expect that this exponent (3 will be asymptotically a constant for enough large n. This is 
a hypothesis and we call it Finite Range Scaling (FRS). 

If the FRS does work, then the infinite n behavior of the susceptibility apart from the finite 
part can be evaluated as 

logx(n->oo) =J^[logx(n) -logx(n-l)] =t(0(p,v))- ( 8 ) 

Using the exponent we estimate the divergent point of the susceptibility and it gives the critical 
dissipation. That is, the infinite n behavior of the susceptibility is controlled by the zeta 
function CC#), which has a pole singularity at /3 — 1. Therefore, the critical 7] is determined by 
the condition (3(p,i] c ) = 1. This method of using finite range systems to evaluate the infinite 
range divergence is the Finite Range Scaling method we have developed. 



3 Effective Ising model of the double-well quantum me- 
chanics 

We move on to the double well quantum mechanics, 

V(x) = -\x 2 + Ax 4 . (9) 

Note that the number of states on a site is infinite in quantum mechanics and T matrix is 
replaced by a bi-local function defined by a bi-local potential W as, 

T(x, y) = e - w(x > y) = (x | U | y) . (10) 
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The operator U is the imaginary time evolution operator and this bi-local function is nothing 
but the Feynman path integral kernel. 

Adding the non-local interactions generated by the path integration of environmental de- 
grees of freedom, we have to compose n sized blocks to make the system 'nearest-neighbor'. 
Then the inter-block T matrix is a multi variable function with 2n variables as follows, 

T = e -^(*i-*». «*•••»») = ( XlX2 ■■■x n \U\ y W2 ■ ■ ■ y n ) . (11) 

Here we introduce a new complete set of states | a n ) and express T matrix with this new base 
system, 

T = e-^-^ b ^ = (a a ...a n \U\b 1 ---b n ) 

dxi ■ ■ ■ dx n dy 1 ■ ■ ■ dy n {ai ■ ■ ■ a n \ x x ■ ■ ■ x n ) (xi ■ ■ ■ x n \ U \ y\ ■ ■ ■ y n ) 
x (yi ■ ■ • y n I h ■ ■ • b n ) 
dx x ■ ■ ■ dx n d Vl ■ ■ ■ dy n r ai (xi) ■■ ■ C (x n )e- w ^- x ^-^ 

*ipbi(yi)---Tpb n (yn), (12) 

where ip a {x) is the wave function of state \a). Our site variables x, y are converted into other 
complete set suffix a, b. 

Now we approximate the states on a site with only two states, that is, the summation 
in a, 6's are just two-fold. This defines an effective Ising model of one-dimension which has 
general interactions between spins, and we can write the multi-local potential as W(er a , cr 6 ) 
where cr a ,cr b = {-1,+1}. 

The simplest version of two-state approximation is to take the completely localized states 
at the left and right well bottoms, where the wave functions are the 5-function. Then the 
converted W(a a , cr h ) function is directly related to the original W(x, y) function as follows, 

W{er a ,(T b ) = W{vcr a ,vcr b ) , v = ^= , (13) 



where ±t> is the position of the double well bottoms. The effective spin interactions in this case 
are given by, 

e2 ~ P £rx ? ? + £ ? akak+l ' (14) 

which is exactly the same form of the normal long range Ising model defined in Eq.([3]) except 
for the total rescaling of the coupling constant rj and the nearest-neighbor interactions due to 
the kinetic term of the quantum system. Then we can easily estimate the critical dissipation 
by referring to the Ising model result calculated in Ref . [HI [T3] as follows: 

rj c = 47rA x r] c [Ising] , (15) 

where we assume that the criticality does not strongly depend on the size of the nearest neighbor 
term. This relation holds for any p and in the Ohmic case (p = 2), for example, we have 

r] c ~ 4vrA x 0.66 . (16) 

This result is quite similar in the A dependence to that obtained by the dilute instanton 
calculation [TB"]. 

Vc = 2n\, (17) 
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and we note that it is larger than this instanton approximation. 

However this selection of states makes difficulty that there appears divergences due to the 
vanishing time discretization slice e. For example the correlation length £ of the system diverges 
as 

^2 eXp (2Al)- (18) 
because of the diverging nearest neighbor interactions. This is not consistent with the fact that 
the quantum mechanics does not need any renormalization of bare parameters. 

Instead, we take linear combination of two states, the ground state |0) and the 1st excited 
state |1) of the double well potential without dissipation (77 = 0), -?/>t,4. = (|0) ± |l))/\/2, which 
are regarded as left or right states corresponding to up or down of Ising spin. We call this 
type of two-state approximation as the ground state approximation. In this approximation the 
correlation length of the system is evaluated as 



log tanh 



1, , /e<T 
■- log tanh I — 



n -1 



(19) 



where 5 is the energy gap between the ground and the 1st excited states. At the vanishing time 
slice, it converges as 

(20) 



lim£(e) 



1 

5 ' 



which is the correct value of the continuum system. 

Hereafter, our aim is to evaluate the plausibility of our BDRG and FRS methods in the 
effective Ising model with the ground state approximation. Note that the effective Ising models 
in this approximation have all possible spin interactions including multi-spin products, which 
might be far away from the normal long range Ising interactions in Eq.([3]). Therefore it is 
non-trivial that FRS method works with such general type of long range Ising models. 

For example, in n = 2, initial T matrix is calculated as follows, 



T = e~ w ^ b ^ = J dx l dx 2 dy l dy 2 ^(^1)^(^2)^(^)^(2/2 
m{x\ — x 2 ) 2 m{x 2 — yi) 2 m(yi — y 2 ) 2 



x exp 



4e 



2e 



4e 



--(V( Xl ) + V(x 2 ) + V( yi ) + V(y 2 )) 



2tt 



~(xx - x 2 ) 2 + (x 2 - yx) 2 + -(yi - y 2 f 



+ ^( { ^ r, - >nf + (x 2 - y 2 ' 2 



(21) 



where e is the discretization step for the imaginary time and a ,b are Ising variables {f ,\.}. 
The first group of terms are the original kinetic terms, the second are the potential terms, and 
the last are dissipation terms given by the CL non-local interactions. 

In the ground state approximation, generally, 2 2n integrations of 2n-dimensions are necessary 
to get the initial T matrix of BDRG. To evaluate these large dimensional integration, we adopt 
the Monte Carlo integration method. We set up random numbers obeying the probability 
distribution function defined by all local terms in the integral, namely, wave functions and 
potential terms. Wave functions are not positive semi-definite and we introduce sign functions 
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to take account of the negative region of the wave functions. Kinetic terms and non-local 
interaction terms are evaluated by random numbers defined above. 

After obtaining the effective spin interactions, we modify it so that the T matrix respects 
the Parity symmetry and mirror symmetry. For example, in n = 2, these symmetries are, 

ty(aia 2 ,6162) = W{a x a 2 , bib 2 ), W{b 2 b x , a 2 a 1 ) = W{a 1 a 2 , b x b 2 ) , (22) 

where \ =\, , \, =t • This symmetrization is important and necessary, since if not, the renormal- 
ization transformation enhances the symmetry breaking components and the resultant renor- 
malized T matrix will be out of reality. We checked the validity and precision of our Monte Carlo 
integration by comparing our results with those obtained by Simpson's integration method 
which can be done for very low n cases. 



4 Results 

We show our numerical results in order. Figj2] shows a typical behavior of the Finite Range 
Scaling exponent f3 versus 77 for n — 5,6, 7. Here we take A = 0.04, the damping rate p = 1.99, 
the discretization step e = 0.9, the number of configurations for Monte Carlo integration =1.28 
million. We calculated with 16 different set of random numbers and evaluated the standard 
deviation of the data which are plotted in the figure as la bar. This initial evaluation is called 
rough estimate. 
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Figure 2: Typical behavior of j3. 



The global behavior of f3 is very smooth and natural. It decreases monotonically with 
respect to 77 and crosses the critical value (3 = 1. According to the FRS hypothesis, the crossing 
point is the critical coupling constant rj c . Note that we need estimates of /3 with enough large 
n, and n = 6 or n = 7 are actually 'large' since it seems they are already stable against n. 
This exhibits a surprising feature that with relatively short range system we can evaluate the 
criticality which can emerge only for infinite range interactions. Such features of FRS method 
is common to the normal long range Ising models. 
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Thus we can obtain the critical coupling constant by FRS method. This holds for all cases 
we calculated with various system parameters X,p,e. 
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Figure 3: Detailed behavior of 0. 
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Figure 4: lcr region of 0. 



After checking the critical region in the above 'rough' data, we increase the Monte Carlo 
data points to be 10.24 million and take 64 different set of random numbers, which we call 
'detailed' data. In fig]3l the results of 'detailed' are plotted with lcr errors, which appears 
consistent with the 'rough' data whose lex region are expressed by curved lines. 

Using the lcr error for each 77 point, we get lcr region of estimate as in figJU Then the 
crossing period of = 1 line with this region gives our estimate of rj c . 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 





A. =0.02 
A =0.04 
A =0.06 
A. =0.10 


i — i — i 

' X ' 


[ 


1 








■ — * — ' 


> 










■ — B 1 










[ 


: 
t 




> 


< 

t 














X 


> 


< 






f 
















- 



11, 



1.2 
1 

0.8 
0.6 
0.4 
0.2 



-0.2 



A =0.02 
A. =0.04 
A =0.06 
A=0.10 


i — I — i 

i X 1 










i — * — i 










' B 1 








* 








> 


< 






X 


) 


< 




) 


< | > 


< * 














i 




- 















0.2 



0.4 



0.6 0.8 



1.2 



0.2 



0.4 



0.6 



0.8 



1.2 



Figure 5: e dependence of 
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Figure 6: e dependence of 
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Next we check the e dependence of T] c . The discretization skip e should be small enough. 
However, the situation is not so simple and we cannot take the vanishing e limit straightfor- 
wardly. The small e means the length scale of the distance between spins becomes small, and 
the long range interaction nature of the system is being lost with fixed finite n. Actually we 
can calculate only rather short range interaction cases around n = 6, 7 because of the computer 
resource limit. Then the too-small e finally breaks the FRS framework. 

There must be an optimized finite e which gives the best value for the physical results 
within the limited calculational resources. Investigating the e dependence of i] c and size of 



their errors in figl5] (p = 1.99, 



n 



6) and figE] (p = 1.99, 



n 



7), we adopt e = 0.9(n = 6) 
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and e = 0.7(n = 7) for major estimates with whole range of parameter values, mainly taking 
account of total stability of calculations. 
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Then we plot the p dependence of r\ c in figJT] (e = 0.9, n — 7). Larger p is equivalent to 
weaker interactions among spins. Therefore ry c should increase according to p. The dependence 
resembles to the case of the normal long range Ising models [UJ. 

However we observe still finite r] c in case of large p > 2, which is quite different from the 
normal long range Ising models defined in Eq.([3]) where r] c is proved to be infinite for p > 2. 
This issue of finite critical coupling constants seen for p > 2 is an interesting problem indicating 
some new features of the model or the important defects of our approximated model. Actually 
our model here does not correspond to a usual simple long range Ising model. Actually we may 
define and evaluate effective Ising spin interactions given by the ground state approximation. 
They are quite different from the simple-minded interactions which are actually realized in 
case of the delta-function states. There appear strong multi spin (a 4 , a 6 etc.) interactions and 
enhanced long range interactions. 

In this sense it is non-trivial that our approximation of the quantum double-well system 
may show up the FRS behaviors just as in the normal long range Ising model, and even more 
the quantum double-well system with long range interactions has the same boundary structure 
at p = 2 as the normal long range Ising model. To conclude something more definite about 
this issue we need to proceed beyond the two-state approximation where FRS will be examined 
with the full system. 

To compare our results with earlier researches, we plot r] c versus A in FigsJS] and M Here we 
take p = 1.99, which is for the future comparison with normal long range Ising calculation with 
FRS since p = 2 there is hard to evaluate due to the boundary Kosterlitz-Thouless transition. 
As for other parameters, we adopt two sets: e = 0.9, n = 6 (FRS(6) in figures) and e = 0.7, n = 7 
(FRS (7) in figures). We show data with error bars which show statistical errors because of the 
Monte Carlo integration part in our method. We see some variation of results between these 
two sets of parameters, which indicate size of systematic errors of our method of calculation in 
this stage. 
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Figure 8: Critical dissipation compared with other results. 



Our FRS results and those of sophisticated Monte Carlo simulation (MC in figures) [T7] 
support complementary parameter regions. At the larger A(> 0.05) area, our results by 2-state 
approximation seems to deviate from the MC, and we consider that it reveals the insuflicientness 
of our approximation of picking up only 2-state at each site in this region. The reason is that 
at larger A region, the well becomes shallower and the energy gaps between the 1st and the 2nd 
excited states are not so large as that between the ground and the 1st excited states. 

On the other hand, the dilute gas approximated instanton results (Instanton in figures) [16] 
look very much consistent with our results at small A(< 0.02) region where both 2-state approx- 
imation and the dilute gas instanton are expected to work well there. The dilute gas instanton 
calculation keeps only the position of instantons and anti-instantons for configurations and 
neglect energy between instantons, that is, non-interacting gas. Then the dissipation effects 
are evaluated for instanton configurations. This approach can be regarded as an another type 
of two-state approximation, and therefore its predictable region resembles to ours. Thus, the 
dilute gas instanton at larger A(> 0.04) is not a good approximation. 

We plot the 2-state approximation results with the delta function wave function (5-wf in 
figures) explained in Eq. (fl6|) though it suffers a serious problem for vanishing e limit. It is just 
above the instanton results and seems better than that. 

Accordingly our method is a good candidate to give the critical coupling constants in the 
intermediate region between the MC plausible region and the dilute gas instanton region. How- 
ever, it seems that within 2-state approximation, our method is not good enough to smoothly 
connect the MC and instanton. We have to proceed to multi-state or full-state calculations, to 
which we will apply the Finite Range Scaling. 

We also plot a non-perturbative renormalization group approach (NPRG-WH in figures) [T^] 
with the local potential approximation, which gives results at very large A(> 0.1) region. This 
approach uses the Wegner-Houghton equation where the effective Wilsonian action consists of 
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Figure 9: Critical dissipation of small A region. 



the fixed kinetic term and the general potential term ignoring all other derivative interactions. 
It is known to give very good results of the susceptibility for large A(> 0.1) caseJ^D]- Ignoring 
the derivative terms in the Wilsonian effective potential means it is not a good approximation 
for configurations with a sharp change of x(t), and it fails for small A(> 0.1) region. In total, 
it looks consistent with other approaches, though there is no other results in the large A(> 0.1) 
region ever. 

In the appendix we list our results of r] c for various parameter values of A,p, e, n. 
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Appendix 

@We list the critical dissipation for various model parameters. The la errors are put in the 
parentheses. 

Critical dissipation r/ c 
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0.63(0.036) 


1.95 


0.049(0.017) 


0.39(0.012) 


0.62(0.018) 


0.79(0.043) 


1.99 


0.060(0.022) 


0.47(0.016) 


0.77(0.023) 


0.95(0.055) 




A ( e = 0.7,n = 6) 


P 


0.02 


0.04 


0.06 


0.10 


1.99 


0.18(0.047) 


0.65(0.026) 


1.12(0.045) 


1.23(0.060) 




A (e = 0.7,n = 7) 


P 


0.02 


0.04 


0.06 


0.10 


1.99 


0.011(0.031) 


0.4(0.024) 


0.80(0.030) 


0.99(0.065) 
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